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Abstract 

In this paper we investigate the electronic and magnetic properties of K A Fe2- v Se2 materials at different 
band fillings utilizing the multi-orbital Kotliar-Ruckenstein's slave-boson mean field approach. We find that 
at three-quarter filling, corresponding to KFe2Se2, the ground state is a paramagnetic bad metal. Through 
band renormalization analysis and comparison with the angle-resolved photoemission spectra data, we iden- 
tify that KFe2Se2 is also an intermediate correlated system, similar to iron-pnictide systems. At two-third 
filling, corresponding to the Fe 2+ -based systems, the ground state is a striped antiferromagnetic (SAFM) 
metal with spin density wave gap partially opened near the Fermi level. In comparison, at half filling case, 
corresponding to the Fe 3+ -based compounds, besides SAFM, a Neel antiferromagnetic metallic ground 
state without orbital ordering is observed in the intermediate correlation range, and an orbital selective Mott 
phase (OSMP) accompanied with an intermediate-spin to high-spin transition is also found. These results 
demonstrate that the band filling and correlation control the electronic state, Fermi surface topology and 
magnetism in K A Fe2- v Se2. 

PACS numbers: 74.70.Xa,71.27.+a,71.10.-w 
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I. INTRODUCTION 



Recently a new iron selenide superconductor K Y Fe 2 _ v Se 2 with T c above 30 K [1] has attracted 
considerable attention for its unique high Neel transition temperature and insulating properties, as 
well as the presence of intrinsic Fe- vacancy ordering [|2j45j], quite different from the other iron- 
based superconducting materials. These unusual properties arise an assumption that K x Fe 2 _ySe 2 
may be a strongly correlated system. KFe 2 Se2 is isostructural to the 122 system, e.g. BaFe 2 As2, 
but chemically is close to FeSe. On average there are 6.5 electrons with equal ratio of Fe 2+ and Fe + 
in KFe 2 Se 2 , rather than 6 ones in the iron-pnictide parent compounds with solely Fe 2+ . Therefore it 
can be regarded as an electron overdoped 1 1 system, in contrast to the underdoped KFe 2 As 2 which 
has 5.5 electrons with equal ratio of Fe 2+ and Fe 3+ , and FeAs which has 5 ones with Fe 3+ . As a 
consequence, only the electron Fermi surface (FS) pockets exist around the M points in KFe 2 Se 2 , 
as observed in recent angle-resolved photoemission spectroscopy (ARPES) experiments [(], 7] and 



in electronic structure calculations [8|,|90- Thus, the FS nesting between the hole pocket around the 
7" point and the electron pocket around the M point, which widely exists in FeAs-based materials, 
is absent in KFe 2 Se 2 compound. 

Due to its unique electronic structure properties, the magnetic properties of K r Fe 2 _ y Se 2 are fo- 
cused on. Because of the difficulty of single crystal preparation for pure AFe 2 Se 2 (A=K, Tl, Rb, or 
Cs), its magnetic properties are mainly studied theoretically, but still a debating issue, as addressed 
in what follows: the local density approximation (LDA) calculations suggested that KFe 2 Se 2 is a 
striped antiferromagnetic (SAFM) order, same to the 1111 and 122 phases of the FeAs-based ma- 
terials Ja]; while some others [8] thought it to be bi-collinear AFM with (n/2, n/2) wave- vector, 
similar to FeTe [|1C 1. However, the LDA calculations suggested that TlFe 2 Se 2 is a checkerboard 
AFM (n, n) order 111 in . and the dynamical spin susceptibility obtained within random phase ap- 
proximation (RPA) also suggested the (n, n) instability in KFe 2 Se 2 11121 . [130. These discrepant 
results show that further investigations on the magnetism are deserved to understand the unique 
properties in AFe 2 Se 2 (A=K, Tl, Rb, or Cs) compounds. 



Many recent experiments Ill4l - l20ll reported the wide existence of the phase separation in 
K v Fe 2 _ y Se 2 materials. A scanning tunneling spectroscopy (STS) experiment demonstrated the 
phase- separated component KFe 2 Se 2 is the parent phase contributing to the superconductivity, 
and Ko. 8 Fei. 6 Se 2 component is an Fe-vacancy order insulator 111411 . implying that pure KFe 2 Se 2 
in the normal state is more possibly a paramagnetic (PM) phase. One may notice that different 



K contents in K r Fe 2 - y Se 2 lead to different band fillings of Fe 3 d orbitals, hence to quite different 
electronic and magnetic properties. Further Chen et al. reported that the electronic states and 
magnetic phase diagrams of K A .Fe 2 - v Se 2 system are closely connected with the Fe valences II 1511 . 
We also notice that various theoretical magnetic configurations obtained within the first-principles 
ations for AFe 2 Se 2 (A=K, Tl, Rb, or Cs) do not include the Coulomb correlation correction 



calcu 

[Si Da lilt!, which implies a weak electronic correlation in AFe 2 Se 2 , in contrast with the intermedi- 
ate electronic correlation in the FeAs-based compounds [21]. Thus a few of questions are urgent 
to be answered: what is the realistic electron filling ? Moreover how does the band filling affect 
the electronic properties and magnetism in K jr Fe 2 _ y Se 2 ? Considering that KFe 2 Se 2 is a possible 
parent phase of superconducting state, we will focus on K l Fe 2 _ y Se 2 system in the absence of 
ordered Fe vacancy at various electron fillings throughout this paper. 

In this paper, to uncover the role of electronic correlation on the groundstate properties of 



AFe 2 Se 2 , we adopt Kotliar-Ruckenstein's slave boson (KRSB) mean field approach [[21 



•220 to 



study the magnetic and electronic properties at different band fillings in K-doped iron selenides. 
Based on our previous LDA calculation results II 1211 . we first present an effective three-orbital 
model for KFe 2 Se 2 , and then determine the ground states of this model at different electron fillings. 
We show that the ground state of K Y Fe 2 _ y Se 2 at fillings of 3/4, 2/3 and a half is a PM metallic phase, 
a SAFM with orbital ordering, and a Neel AFM one without orbital ordering in the intermediate 
and strong correlation regimes in addition to an orbital selective Mott phase (OSMP) related with 
an intermediate-spin to high-spin transition, respectively, showing that the band filling controls 
not only the FS topology, but also the electronic structure and magnetic properties of Fe-based 
superconducting materials. The rest of this paper is organized as follows: a three-orbital tight- 
binding model and the multi-orbital slave-boson mean-field approach are presented in Sec. II; the 
numerical results and discussions are shown in Sec. Ill; the last section is devoted to the remarks 
and summary. 



H. THREE-ORBITAL TIGHT-BINDING MODEL AND SLAVE BOSON APPROACH 



Based on our previous electronic structure calculations [12], we find that FS is mainly con- 
tributed by three h„ orbitals, thus the system can be described by a three-orbital model, similar to 
iron pnictides ll24M27ll . We extract an effective three-orbital tight-binding model from our LDA 
band structures. The tight-binding model Hamiltonian for the three-orbital model in the momen- 



turn space is described as, 



(1) 



kacr 



where T a/3 (k) is the kinetic energy term, e a denotes the on-site energy of the a orbital, and /j. is the 
chemical potential determined by the electron filling. The three-orbital tight-binding fitting of the 
Fe-3J bands is displayed in the solid lines, in comparison with the original five bands J12Q in the 
dot lines, as shown in Fig. 1. It is obviously found that the band structures in KFe 2 Se 2 are similar 
to that of LaFeAsO |28|], only the position of Fermi energy E F is shifted. Therefore this model can 
describe both the FeSe-based and FeAs-based systems through changing the chemical potential. 
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FIG. 1. (Color online) Left panel: the band structures of the Fe-3d orbitals obtained by the LDA and its 
three-orbital tight-binding fitting. Dot lines are the LDA five bands and solid lines the fitted effective 
three bands. Right panel: the corresponding Fermi surface for parent material KFe2Se2 at the first Brillouin 
zone obtained by the three-orbital tight-binding model. 



The tight-binding parameters of the three-orbital model for KFe 2 Se2 are listed in the follow- 
ing. The on-site energies measured from the Fermi energy for the three orbitals are (e\, 62, £3)= 
(-511.92, -511.92, -341.63), respectively, in units of meV. Here orbital indices (1, 2, 3) indicate 



the d xz , d yz , and d xy components, respectively. Similar to BaFe 2 As2 fl29J], the orbital-dependent 



kinetic energy matrix elements are expressed in terms of the inter-obrital and intra-orbital hopping 



integrals as follows, 

T 11/22 = 2t l xjy cosk x + 2t l y ) x cosk y + 4t^cosk x cosk y 

±2t x \{cos2k x - cos2k y ) + 4t xxy/xyy cos2k x cosk y 
+4t xyy/xx) ,cos2k y cosk x + At l xxyy cos2k x cos2k y , 

T 33 = 2t 33 (cosk x + cosk y ) + At' 3 cosk x cosk y + 2t 33 x (cos2k x + cos2k y ) 
+4t 3 xxy (cos2k x cosk y + cos2k y cosk x ) + At 33 xyy cos2k x cos2k y , 

T n = At l2 sink x sink y + 4t xxy (sin2k x sink y + sin2k y sink x ) 
+At xxyy sin2k x sin2k y , 

T l3/23 = ±2it 13 sink x/y + Ait X3 cosk y i x sink xjy 
+4itl xy co sk y / x sin2k x / y , 

The intra-orbital and inter-orbital hopping parameters up to the fifth nearest-neighbor for the fitting 
of the three-band structure in Fig. 1 are shown in Table /. 

TABLE I. The nonzero matrix elements of intra-orbital tf a and inter-orbital t°f hopping parameters up to 
fifth neighbors of the three-orbital tight-binding model through fitting the band structures. All the parame- 
ters are in units of meV. 
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Considering the Coulomb interaction, in addition to the kinetic term in Eq. (1), we describe the 
electronic interaction part of the multi-orbital Hamiltonian as follows, 

Hi = n ia \n iai + U' ^ n i^ n ipi + W - J H ) ^ 



i,a i i,cr 



i i 
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where U(U ) denotes the intra-(inter-)orbital Coulomb repulsion interaction and 7 H the Hund's 
rule coupling. Considering the rotation symmetry of the system, we adopt f/' = C/-27 H - 

We notice that in FeAs-based materials, the electronic filling for the present three-orbital model 



is only at two-third filling, i.e. n=4, the correct FS can be reproduced 11261 - 12 811 . In KFe 2 Se 2 
compound, the Fe ions have two kinds of valence, Fe 2+ and Fe + with equal ratio, and the average 
electron number is 6.5 per site, different from Fe 2+ with 6 electrons in FeAs-based parent materials. 
Consequently, pure KFe 2 Se 2 should be at a filling of three quarters (i.e. n=4.5) at which the FS 
can be reproduced correctly in the present three-orbital model. The obtained FS of KFe 2 Se 2 is 
plotted at k z =0 in Fig. 1, with four electron-like FS pockets. The hole-like FS pockets at 7" are 
absent. Such a FS topology is in agreement with the recent ARPES experiments fls |7|] and the 
band structure calculations |8, sj]. To explore the band filling dependence of the electronic state, 
we also extend the band filling to 2/3 and half, which corresponds to the Fe 2+ and Fe 3+ systems, 
respectively. 

KRSB mean field approach is known as one of the effective methods to treat the wide elec- 
tronic correlation in the many-body systems. Here, we extend the single-orbital KRSB mean-field 



method to the multi-orbital situation 1121 . 



2311 and apply it on K v Fe 2 _ y Se 2 . We introduce new 
fermion operators^ slaved by boson operators e„ p ia(T , d iaa - p(Ty , b ia , t ia/3(T , r i(Ta(Tf>0 - y , q ia , u ia(Tp(Ty , 
Viaa- and s, which represent the empty, single occupation with a orbital and <x spin, double occupa- 
tion with [3 orbital spin <jp and y orbital spin cr y , double occupation with two electrons in orbital 
or, triplicate occupation with two electrons in orbital a and one electron in orbital with spin <x, 
triplicate occupation with each electron in orbital a, f3 and y with spin cr a , cr^, and cr y , quaternity 
occupation with two electrons in orbital /? and y, quaternity occupation with two electrons in or- 
bital a, one electron in orbital fi with spin erg and one electron in orbital y with spin cr y , fivefold 
occupation with four electrons in orbitals and y, and one electron in orbital a with spin cr, and 
sixfold occupation, respectively. 

The completeness of these boson fields gives rise to the normalization condition as follows: 

e l e ' + Pjgo-Piacr + b\cPia + / j ^lacrpCTy^'^^y + AcTaCr^^a^^y 

+ Z t tfr t **°- + Z q te qU * + Z U l^ 7 U i^y + J] V i<r V ™ + A S > = 1 ( 3 ) 
a,fi,cr a a,crp,cTy a,<r 

Projecting the original Hamiltonian Eqs. (1) and (2) into the new slave-boson representation, the 



6 



multi-orbital Hubbard model Hamiltonian is described as: 

<xJ3,y 

i,j,a,p,cr iacr i,a i,a 

+(U' - J H ) £ + £/' J] + (U + 2£/' - /„) ^(^W + 

+/ff ^(^•ggo-^yag- + tiPa<jt\ yaa ) + 3(C/ - Jh) ^ r L a tr jS (r r r icr a crpcr y 
i,a,cr i,cr 

a,p,y ct,p,y 

+(3U' - J H ) J] r la a <T P a/aa a ^fy + 2 ( U + 2U ' ~ J h) J] vial** + Jh Z^W'y + 
i,a,<r i,a i,a 

aJ3,y ct,p,y 

+(U + 5U'- 3J H ) J] "L-^"^ +(U + 5U'- 2J H ) Y u] mTpif u ia(Tp(fy 

i,a,cr i,a,cr 

+(2U + W - AJ H ) J] v] ao .v ia(T + (3U + 12U' - 6J H ) J] s]*i (4) 

i,a,cr i 

where the renormalization factor 

Ziacr = (1 — Qiacr) 2 ^iacrQj a(T (5) 

a,p,y a,/3,y P+cy 

%acr = e\ Pica- + p\ ad b ia + ^ P^diycra' + Plycr'dipcr'cr + ^ b^tip aa 

cr' cr' p 

afi,y ct,p,y ct,p,y 

0~',0~" 0~' Cr, Cr',Cr" 

"^hya&QiP + hpa&QiJ + 2-i iyPcr' U 'y a ~<*°~P + 2_j hpycr' U fiVyVa + 1ia Viao ~ 
cr' cr' 

+ Z U Ua<r/i^' + J] U la y a a V iy<r' + V lcr S i ( 6 ) 
cr' cr' 

The corresponding Fermion number constraint for a orbital with spin <x reads: 

Qiacr ~ f[ acr fiacr (7) 

where 

Gi'acr = Pjao-Piacr + ^fo^ia + ^ Apcr y cr a ^WyO- a + ^ Aycr^^iy^^P + ^ r ]cr a crpcr y r i"-a^y 

cr' cr' cr',cr" 

P*a 



+ ^iaPcrtiaPo- + ^ ^ipajipacr + ^ <7^<?,£ + ^ U ]ycr„cr p U iyo- a crp + ^ U ](3cr y cr a U iP°-y°- a 
P,cr p P cr' cr' 

P*a 



Averaging the boson operators in Eq.(4)-(8), we can obtain an effective mean filed Hamiltonian, 
hence the total ground-state energy. The original fermions are guaranteed by the constraints Eq.©, 
which are implemented by means of the corresponding generalized Lagrange multiplier method. 
In order to determine the stable magnetic ground state, we minimize the total energies for different 
magnetic configurations based on the pattern search method. To simplify the calculations, various 
symmetries should be utilized. For instance, in the AFM situation with two sublattices, the single 
occupation probabilities p^ in sublattice B are identical to p^ in sublattice A, double occupation 
probabilities , =d\^, etc.. 



HI. RESULTS AND DISCUSSIONS 



In this section, we present main numerical results on the electronic and magnetic properties 
within the three-orbital model for iron selenide systems. The fillings of three quarters and two 
thirds, as well as half-filling, are all considered for comparison. 



A. Three-quarter filling case 

We firstly consider the electron filling n=4.5 case, corresponding to pure KFe 2 Se 2 compounds. 
Note that we adopt the hole representation for convenience within KRSB mean field approach 
throughout this paper. Thus the particle number at the three-quarter filling case is 1.5 within 
the three-orbital model. Taking into account several types of magnetic configurations with high 
symmetry, the PM, ferromagnetic (FM), Neel AFM and SAFM cases, we find that only the PM 
phase is the most stable at a filling of three quarters when U increases from up to 5 eV. Because of 
the mixing valence of Fe 2+ and Fe + in KFe 2 Se 2 , the system is not a magnetic ordered and insulating 
state in the homogenous phase. Since the electron fillings of both Fe 2+ and Fe + are away from half 



filling, unlike mixing valent Na 5 CoO 2 [30], KFe 2 Se 2 does not form charge ordering. Previous 



work on charge susceptibility shows that the Coulomb interaction suppresses the charge instability 



3 111 , thus it is hard to form charge ordering in KFe 2 Se 2 . 
Within the present KRSB framework, the dependence of various boson occupancies on the 
Coulomb interaction U at / H =0.25£/ is plotted in Fig. 2. It is clearly found that the empty oc- 
cupation e, single occupation p 3 and double occupation b 3 with orbital xy are dominant for small 
U. With the increase of the Coulomb interaction, the empty occupation e and triple occupation 
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FIG. 2. Dependence of boson occupancies on Coulomb interaction U at an electron filling of n-4.5, and 
/H=0.25t/. Note that in the PM case, do a - ^a-u = ^«1T> d\ a - d a ^=d a n. 
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FIG. 3. (Color online) Orbital occupations (left panel), and renormalization factor (right panel) as a function 
of the Coulomb interaction U at n-4.5, and Jh=0.25U. 

b 3 decrease sharply, while du, d 12 and d ]3 (di a =d a ^=d a n) increase due to the increasing of the 
Hund's rule coupling and the Coulomb interaction. This shows that the system is a PM and the 
magnetic moment of each Fe spin increases with the lift of U and rather than a nonmagnetic 
one. 

The orbital occupations as a function of the Coulomb interaction are given in the left panel 
of Fig. 3. It is obviously found that the electrons are occupied at the lower energy xz- and yz- 
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FIG. 4. Projected densities of states are plotted for U=0 (a), 1.0 eV (b), 3.0 eV (c), and 5.0 eV (d) at n=4.5, 
respectively, with Jn=0.25U. Inset: total densities of states. 



orbitals. And there is no orbital polarization between the xz- and _yz-orbitals. With the increase 
of Coulomb interaction, a fraction of electrons in the xz- and yz-orbitals transfer to higher energy 
xy-orbital since larger J H enhances the effect of the Hund's rule. Thus the high energy xy orbital 
occupies more electrons as J H increases. The renormalization factors of each orbital as a func- 
tion of the Coulomb interaction U are shown in the right panel of Fig. 3. With the increase of 
Coulomb interaction, the bandwidths of the three orbitals become narrower and narrower, and the 
renormalization factors become small. We notice that in the PM phase, the degeneracy of orbital 
xz and yz is not removed, and the renormalization factors of bandwidths in different orbitals are 
nearly the same. When U is 3, 4 and 5 eV, the renormalization factor Z is about 0.85, 0.8 and 
0.76, respectively, giving rise to the band mass of original fermions m b /m = 1/Z 2 ~ 1.38, 1.56 
and 1.73, respectively, and yielding a mass renormaization factor of about 2 at U=5 eV. We ex- 
pect that the disorder effect and the spin fluctuations beyond the KRSB mean field approximation 
will further narrow the bandwidths and enhance the effective mass. This band narrowing fac- 
tor is comparable with the experimental ones from ARPES (32D and de Haas- Van Alphen [33], 
indicating that KFe 2 Se2 lies in the intermediate correlation region, similar to the FeAs-based sys- 
tems. While in the presence of the Fe-vacancy, the measured renormalization factor is about 6.1 
in (Tl,Rb) r Fe 2 - v Se2, greatly larger than 2. This shows that most probably the insulating proper- 
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ties in K v Fe2- v Se2 systems are induced by the ordering of Fe vacancies, rather than by the strong 
electronic correlation. 

The projected densities of states (PDOS) of KFe2Se2 for the Coulomb interaction U =0, 1 , 3 and 
5 eV are shown in Fig. 4. The finite DOS at FS shows a pseudo-gap-like structure, and there is no 
obvious van Hove singularity. Nevertheless in FeAs-based compounds, a considerable van Hove 
singularities in the DOS is attributed to the FS nesting with wave vector Q=(n, 0) 04I1 . Therefore, 
the absence of the van Hove singularities in the DOS of KFe 2 Se 2 may be the consequence of the 
breakdown of the FS nesting, which arises from the fact that there lacks the hole pockets around T 
point. With the increase of the Coulomb interaction from to 5 eV, the bandwidth becomes narrow, 
W ~ 2.5 eV at U=3 eV, smaller than 3.2 eV at U=0, showing that the system is a correlated bad 
metal and is in the intermediate correlation region with W~U. 

Through the analysis above, we find that at 3/4 filling, the K x Fe2-ySe2 system, corresponding 
to KFe 2 Se2, is a PM phase. We also clarify that KFe 2 Se 2 lies in intermediate correlation region, 
similar to the FeAs-based compounds. 

B. Two- third filling case 

On the other hand, the electron filling at n=4 corresponds to Fe 2+ -based compounds, such 
as FeSe and LaFeAsO, etc.. The dependence of boson occupancy probabilities on the Coulomb 
interaction is also shown in Fig. 5. Comparing with n=A.5 case, we find that with the increase 
of U, the system transits from a PM metallic phase to a SAFM metallic one at a critical point 
U c x\.2 eV. In the PM phase, the single occupations in the xy-orbital p 3 and empty occupancy e 
are dominant. Other multiple occupation states also contribute small but finite weights. While in 
the SAFM phase, the single, double and triple occupations with the same spin alignment, p, d and 
r, are dominant. The spin singlet and small spin states contribute very little. With increasing the 
Hund's rule coupling, the empty and single occupations continuously decrease, while the double 
and triple ones considerably increase. This behavior arises from the fact that the increase of 
Coulomb interaction favors the formation of the SAFM phase, and large Hund's rule coupling 
favors large spin state. 

The gap opening behavior of the SDW states in the multi-orbital FeAs-based system is an inter- 
esting but unsolved topic. In order to resolve the behavior of the SDW gap opening in K r Fe2_ > Se2, 
we present the band dispersions of the PM and SDW states in Fig. 6. We find that in the SDW 
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FIG. 5. Dependence of various boson occupancies on the Coulomb interaction U at n-4, and Ju=0.25U. 




FIG. 6. Band dispersions of the PM and SDW states along the high symmetry points, r (0,0), Y (0, n), M' 
(n/2,7t), X' (tt/2,0), in the folded BZ with U=0 and 2 eV, and J B =0.25U at n=4. 

states, only partial SDW gaps open near E F in comparison with PM states, which is the con- 
sequence of the the band narrowing and spin splitting, similar to FeAs-based compounds. The 
partial opening of SDW gap indicates the system is a bad metal. 

The dependence of orbital occupations and magnetic moment of each orbital on the Coulomb 
interaction U is also shown in Fig. 7. With the increase of U, once the system enters the SAFM 
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FIG. 7. (Color online) Dependence of orbital occupations and magnetic moment of each orbital (left panel), 
and renormalization factor (right panel) on the Coulomb interaction U at n-A, and /h=0.25£/. 



phase, there is obvious orbital polarization between xz and yz orbitals together with n xz < n yz in the 
electronic representation. And the magnetic moments on the xz and yz orbitals are different with 



35|]. The renormalization 



m X7 > m yz . These are consistent with the previous LaFeAsO results [|2lL 
factors Z displayed in the right panel of Fig. 7 show that the Z xz ^, Z >z1 - and Z^. T become smaller 
and smaller with Z xz / yz ^< when U increases. This indicates that the SDW gap opening mainly 
occurs in the xz/yz orbitals. 



The evolution of the PDOS with Coulomb interactions U is also displayed in Fig. 8. We find 
that at [7=2.0 eV and 3.0 eV, there exhibit obvious spin polarizations in SAFM phase due to the 
breakdown of spin symmetry. The xz- and yz- orbitals dominate the FS, while the xy orbital mainly 
lies far from the Fermi levels. This also shows that the xz/yz orbitals determine the electronic 
properties near E F , and mainly involve with the formation of the SDW states, consistent with the 
result of renormalization factor in Fig. 7. 



Since the electronic correlation is intermediate, one expects that U is larger than 2 eV in iron 
pnictides and selenides; thus at 2/3 filling, the SAFM state with orbital ordering is stable over a 
wide electron correlation range, from intermediate to strong correlations, which is similar to the 
FeAs-based systems. 
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FIG. 8. Projected densities of states are plotted for U=0 (a), 1.0 eV (b), 2.0 eV (c), and 3.0 eV (d) at n-A, 
respectively, with 7h-0.25£/. Inset: total densities of states. 

C. Half-filling case 

When turning to the half-filling case, which corresponds to the Fe 3+ -based systems, such as 
KFeSe 2 and FeAs compounds, etc., we find that the phase diagram becomes much richer. It is 
found that there exist three critical points when U increases: the system transits from a PM metal 
to an SAFM metal at U CI , from an SAFM metal to a Neel AFM metallic phase at U C2 , and from an 
AFM metal with intermediate-spin state to an AFM OSMP with high-spin state at U C3 , as shown 
in Fig. 9. These phases will be addressed in detail in the following. In comparison with n=4.5 
and 4 cases, we find that besides the PM and SAFM phases, the Neel AFM metallic phase appears 
in slightly large U region at C/ C2 ~0.8 eV. Only when in the narrow Coulomb interaction region 
at t/ Cl ~0.65<f/< U C2 , the SAFM phase is stable, as the dependence of the boson occupancies on 
Coulomb interaction U showing in Fig. 9. With the increase of the Coulomb interaction, the triple 
occupation sharply increases, while the other high occupations are relatively small. This shows 
that the system undergoes a spin state transition with the increase of the Coulomb interaction and 
the Hund's rule coupling, as we see the low-spin (5 = 1/2) to intermediate- spin (5 = 1) transition at 
U C2 in Fig. 10. 

On the other hand, the dependence of the orbital occupations and the magnetic moments of 
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FIG. 9. Dependence of boson occupancies on Coulomb interaction U at half-filling n=3 and Jn=0.25U. 
Three dashed vertical lines denote the phase boundaries of the PM-SAFM, SAFM-Neel AFM, and Neel 
AFM-OSMP transitions at U Cl , U C2 and U C3 , respectively. 
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FIG. 10. (Color online) Orbital occupations and magnetic moment of each orbital (left panel), and renor- 
malization factor (right panel) as a function of Coulomb interaction U at half-filling n=3, and Ju=0.25U. 
The dashed lines indicate the phase boundaries U Cl , U C2 and U C3 

three orbitals on Coulomb interaction U is plotted in the left panel of Fig. 10. We find that in the 
SAFM phase, a small orbital polarization appears with n yz > n xz in the present hole representation, 
i.e. n xz > n yz in the electron representation. This supports the itinerant orbital ordering in the parent 
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phases of iron pnictides [23, 35]. Meanwhile the magnetic moments on different orbitals possess 
m yz < m xz with total magnetic moment m tot <lfi B - The system lies in a low-spin state. However, 
in the Neel AFM metallic phase, there is no orbital polarization, since in the presence of the 
spin-orbital coupling, the preserved spin rotational symmetry does not lift the orbital degeneracy 
or break the orbital symmetry. When [7>[7 C3 , the system enters the OSMP phase, as displayed 
in the left panel of Fig. 10. All of the three orbital occupations are nearly equal to 1. The 
magnetic moment per orbital steeply increases, and the total magnetic moment m tot is larger than 
2fi B - Meanwhile, the system undergoes an intermediate- spin (5 = 1) to high spin (5 =3/2) transition 
with the increasing of [7, and thus enters a high-spin state when U>U C3 . 

The right panel of Fig. 10 displays the renormalization factor of each orbital as a function of 
Coulomb interaction [7, which is in sharp contrast with the n=4.5 and 4 cases. In the PM and 
SAFM phases, all the renormalization factors Z xz , Z yz and (Zx^kZ^) smoothly decrease with 
increasing [7, indicating the bandwidths become narrow due to the increase of the Coulomb cor- 
relation. When U>U C2 , the system enters the Neel AFM state. In this situation Z^, gradually 
decreases, while Z xz and Z yz considerably change with [7, suggesting that the variations of orbital 
and magnetic states mainly occur in these two orbitals. The lift of Z xz and Z yz implies the band- 
widths of the two orbitals anomalously increase, which is attributed to the fact that the exchange 
splitting of spin-up and spin-down subbands of the xz- and yz-orbitals increases with the increase 
of [7, leading to the total bandwidth broadening, which is also seen in Fig. 1 1 in what follows. 

We present the PDOS of four typical phases, including PM, SAFM, Neel AFM metal and AFM 
OSMP, in Fig. 1 1 for Coulomb interaction [7=0, 0.75, 1.0 and 1.5 eV, respectively. It is found that 
in relatively small [7 region, i.e. U< U C2 , the xz/yz orbitals dominate the FS in both the PM and 
SAFM phases. Different from Fig. 1 1(a), in the SAFM phase with £7=0.75 eV in Fig. 1 1(b), the 
orbital degeneracy between xz and yz orbitals is lifted, consistent with the orbital polarization in 
the left panel of Fig. 10. Meanwhile in the Neel AFM metallic state with [7=1.0 eV in Fig. 1 1(c), 
the weight of the xz and yz orbitals in the FS is greatly suppressed; but the spin splitting becomes 
large, consistent with the intermediate spin configuration in Fig. 10. Interestingly, we find that 
when U>U C3 ~ 1.1 eV, an orbital-selective Mott transition occurs, i.e. an OSMP emerges. In this 
situation, only the broad xy orbital contributes to the FS, while the narrow xz and yz orbitals are 
insulating and sank below the FS, as seen the PDOS with [7=1 eV in Fig. 1 1(d). 

Therefore, at a half filling, the system is mainly a Neel AFM state without orbital ordering in the 
intermediate and strong correlations, in addition to an OSMP phase related with an intermediate- 
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FIG. 11. Projected densities of states are plotted for U=0 (a), 0.75 eV (b), 1.0 eV (c), and 1.5 eV (d) at 
half-filling n=3, respectively. Jn=0.25U. Inset: total densities of states. 

spin to high-spin transition. Compared with the other band fillings, the different magnetic phase 
diagrams suggest a band-filling controlling magnetism scenario in the iron selenide systems. 

IV. DISCUSSION AND SUMMARY 

The scenario of the orbital selective Mott transition in K A Fe 2 _ySe2 could be understood in a 
sketch of OSMP in the presence of magnetism shown in Fig. 12 through tuning the band filling and 
the correlation strength. When tuning to the half-filling, with the increase of Coulomb interaction 
and Hund's rule coupling, the spin-up and spin-down subbands split with each other, and the spin- 
up xl- and yz-orbitals are filled, sank below E F and become insulating; while the xy-orbital remains 
across E F and is conducting, giving rise to an AFM OSMP. We notice that in comparison with the 
bandwidths of K,.Fe 2 _ySe 2 , the U C3 for the occurrence of the OSMP is small. This may arise from 
the following two reasons: (1) in the parent phase of Fe-based superconductors, the Hund's rule 
coupling is large with J H =0.25U. If small values 7#=0.15t/ or 0.1U are adopted, U C3 will reach 
a relatively large values up to 2 or even 3 eV; (2) U C3 is obtained for the OSMP with AFM order 
in the present study, it will become larger in the PM situation. Since KFeSe 2 possesses not only 
tetrahedra crystal field but also large magnetic moment, it may be a potential candidate for such an 



17 



OSMP phase. We expect that further increasing the Coulomb interaction U will drive the system 
to a magnetic insulating state accompanied with a metal-insulator transition. 
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FIG. 12. Sketch of sublattice band structures in the orbital-selective Mott phase with Neel antiferromag- 
netism for half-filled three-orbital model in K A Fe 2 _ v Se 2 . 



We notice that different from our PM results, the first-principles electronic structure calcula- 
tions suggested the ground state in KFe 2 Se2 is the SAFM ordering [9], or the bi-collinear AFM or- 
dering resulting from the interplay among the nearest, the next-nearest, and the next-next-nearest- 
neighbor superexchange interactions mediated by Se 4p-orbitals in _AFe 2 Se 2 (A=K, Tl, Rb, or Cs) 
| 8J1. In addition, some other iron-based materials, such as FeSe I 36L 
[38 -Ek 



370, LiFeAs |38|], and KFe 2 As 2 



0], etc., have no magnetism at all as observed in the experiments, but AFM is obtained in the 
LDA calculations. The discrepancies among these materials suggest the magnetism is sensitive 
to the electronic properties or lattice distortion Q38D . The LDA methods usually so overestimate 
the magnetic moment and AFM ordering but omit some spin fluctuations of the system due to the 
intermediate electronic correlation that a proper treatment on the electronic correlation in these 
FeSe-based compounds may be important for understanding its magnetic ground state. In our 
study, we deal with the electronic correlation within the framework of the KRSB approach which 
is verified as an effective approach to treat the electronic correlation ranging from weak through 
intermediate to strong ones. Our PM groundstate result, not non-magnetic one, indicates that there 
exists strong spin fluctuation which is regarded as the superconducting pairing mechanism, thus 
KFe 2 Se 2 is a potential candidate of the parent phase for superconductor without doping. More- 
over, a recent STS experiment demonstrated the existence of phase separation and distinguished 
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the contributions: the KFe 2 Se 2 component contributes to superconductivity, while the Ko.gFei.gSea 
Fe- vacancy ordering component to the insulating properties [|14j], suggesting that pure KFe 2 Se 2 is 
more possibly a PM phase when in the normal state. These observations are consistent with our 
present results. 

In summary, starting with an effective three-orbital model for newly found KFe 2 Se 2 , we have 
shown that the ground state of KFe 2 Se 2 at three-quarter filling is a paramagnetic metallic phase. 
We also suggest a possible OSMP phase through tuning the band-filling in K v Fe 2 _ y Se 2 . Our results 
demonstrate that the band filling plays a key role in the electronic and magnetic properties of 
K Y Fe 2 - y Se 2 . Since the phase separation widely exists in K x Fe 2 _ v Se 2 materials, future experiments 
are expected to confirm the magnetic ground state of KFe 2 Se 2 . In addition, the influences of the 
band filling and correlation on the electronic structures and magnetic properties in the presence of 
Fe vacancies for K x Fe 2 _ y Se 2 is an interesting topic, thus further theoretical works are expected. 
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